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Abstract. This is a review of spin-dependent (polarization) properties of multicom- 
ponent exciton-polariton condensates in conditions when quasi-equilibrium mean- 
field Gross-Pitaevskii description can be applied. Mainly two-component (spin states 
±1) polariton condensates are addressed, but some properties of four-component 
exciton condensates, having both the bright (spin ±1) and the dark (spin ±2) com- 
ponents, are discussed. Change of polarization state of the condensate and phase 
transitions in applied Zeeman field are described. The properties of fractional vor- 
tices are given, in particular, I present recent results on the warping of the field 
around half-vortices in the presence of longitudinal-transverse splitting of bare po- 
lariton bands, and discuss the geometrical features of warped half-vortices (in the 
framework of the lemon, monstar, and star classification). 



1.1 The Gross-Pitaevskii equation 

Currently, most theoretical descriptions of exciton-polariton condensates ob- 
served [1,2,3,4,5 in incoherently excited semiconductor microcavities are 
based on the Gross-Pitaevskii equation (GPE). When the polarization of the 
condensate is of interest, this equation can be generically written as 

where the order parameter ip(r,t) of the condensate is a complex 2D vector 
function of the 2D position in the microcavity plane r and time t. Alterna- 
tively, one can expand i/> on the circular polarization basis 

, x + iy x-iy 

^ = -^-^ +1 + -^-^ ll (1.2) 



and obtain the coupled GPEs for two circular components ip±% (see Eqs. (|1.32b .b) 
in Section lL4l) . 
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GPE is used in two main flavors, strongly non-equilibrium GPE and 
quasi-equilibrium GPE, that treat the energy relaxation in two extreme ways. 
It is completely neglected in the former, and it is considered to be essential 
in the later. Mathematically, these approaches differ in Hamiltonian density 
H(ip, ip*)'- it is complex in the former and it is real in the later. Each approach 
has its benefits and drawbacks. 

The imaginary part of the Hamiltonian for non-equilibrium GPE [6j[7j|8j|9] 
is given by the difference of income and escape rates of exciton-polaritons into 
and out of the condensate. While the escape rate is given by the reciprocal 
radiative life-time of exciton-polaritons and is independent of the particle den- 
sity, the income rate is a non-linear function of it. The nonlinearity is essential 
to stabilize the solution and it appears due to the depletion of an incoherently 
pumped reservoir. This approach resulted to be quite successful in modeling 
the experimental data on condensate density profiles for spatially nonuni- 
form exciton-polariton condensates. On the other hand, it cannot describe 
the spontaneous formation of linear polarization of the condensate — the fact 
that is quite unfortunate since the observation of spontaneous linear polariza- 
tion is one of the direct experimental evidences of Bose-Einstein condensation 
of exciton-polaritons [UE] . A workaround is to add the Landau-Khalatnikov 
relaxation into (jl.ll) , which is equivalent to adding an imaginary part to time 
t. But this relaxation is rather artificial because it changes the number of 
particles in the condensate. 

In what follows, we consider the opposite limit assuming a fast relaxation 
of exciton-polaritons, so that they reach quasi-equilibrium, even with the tem- 
perature that can be different from the lattice one. The balance of income and 
outcome rates produces some steady-state concentration of exciton-polaritons, 
that can be defined, as usual, by introducing the chemical potential fi. The 
coherent fraction of condensed particles can be described by the traditional 
GPE with real H(ip,ip*), 



Here T is the density of the kinetic energy, 7^i n t describes interaction between 
the particles, H' stands for some possible perturbations, and n = ip* ■ ip is 
the exciton-polariton density. 

The kinetic energy of exciton-polaritons in planar microcavities depends 
on the orientation of vector ip with respect to the direction of motion. Near 
the bottom of lower polariton branch one has 



H = T - M« + "Hint + 



(1.3) 




(1.4a) 



h 2 dip + i dtp 2 h 2 di/j + x dip -i 2 



(1.4b) 



mi dz* dz mt dz* dz 



where mi and nit are the longitudinal and transverse effective masses of po- 
laritons, and the complex derivatives 
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d_ _ 1 ( d_ _ .d_\ d _ 1 / d .d_\ 

dz ~ 2 Vfe ' ~ 2 Vfe ' 

are used. The vector i/> is proportional to the in-plane electric field vector 
of exciton-polariton mode. According to (jl.4l) . the frequency of transverse 
electric (TE) mode with in-plane wave vector k _L %p is hk 2 /2mt, while for the 
transverse magnetic (TM) mode with k\\i/) the frequency is hk 2 /2mi. (The 
same bare frequencies of both modes at k = are removed from (|1.4|) .) 

The polariton-polariton interaction is also anisotropic: it depends on mu- 
tual orientation of i\) and ip* . One can construct two quartic invariants from 
these two vectors and 'Hint is given by 

H int = \(U - EW- iP) 2 + ±Ui\1>* x ip\ 2 (1.6a) 

= |l/o (|V + i| 4 + l^-il 4 ) + (Uo - 2C/ 1 )|^ +1 | 2 |V-i| 2 . (1.6b) 

It is seen that Uq is the amplitude of interaction of polaritons with the same 
circular polarization (with the same spin), and Uq — 2U\ is the amplitude of 
interaction of polaritons with opposite circular polarizations (opposite spins) . 
These quantities are denoted by ot\ and «2 in some papers. The constant Uq 
is positive and can be estimated as ~ Stag, where £b is the exciton binding 
energy and as is the exciton Bohr radius. The interaction of exciton-polaritons 
with opposite spins depends substantially on the electron-electron and hole- 
hole exchange processes and is defined by the electron and hole confinement 
within quantum wells and by the number of quantum wells in the microcavity. 
As a result, the value of Ui is sensitive to the microcavity geometry. 

To end this section it is important to mention the limitations of any GPE 
in application to the condensates of exciton-polaritons in microcavities, or to 
condensates of any other bosonic excitations that have a finite radiative life- 
time. Due to interference of light emitted from different parts of condensate 
there appears dissipative long-range coupling in the system. Most importantly, 
the escape rate becomes dependent on the symmetry of the condensate wave- 
function and this favors the formation of particular long-living many-particle 
states [IU] ■ These effects cannot be properly treated in the framework of Gross- 
Pitaevskii equation (11.11) . 

1.2 Polarization and effects of Zeeman field 

The interaction energy (jl.6p of the polariton condensate is polarization de- 
pendent. While the first term in (|1.6b ) does not depend on polarization and is 
simply proportional to the square of the polariton concentration n = ip), 
the second term in (|1.6a ) is sensitive to the degree of the circular polarization 
of the condensate. For Ui > the interaction energy is minimized when the 
second term in (jl.6a ) is annulated, which is achieved for polarization satisfy- 
ing tp* x ip = 0, i.e., for the linear polarization. On the other hand, in the case 
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Ui < the minimum is reached for the circular polarization of the condensate, 
when ip* x ip = ±in. 

So, there is qualitative change in the ground state of the condensate when 
U\ changes sign [TT] . 

(i) U\ > 0. The ground state is characterized by two angles, the total phase 
angle 9 and the polarization angle 77. These angles are defined from the 
Descartes components of the order parameter ip x = ^Jn e cos r\ and ip y = 
y/ne lS sin rj. The circular components are then ip±i = \Jnj2 e I ^ =F?) - ) . There 
are two broken continuous symmetries and, consequently, the excitation 
spectrum consists of two Bogoliubov branches. The sound velocities for 
these branches at mi — m t — m* are v = -y//i/m* and v\ = y^nUi/m*, 
where /i = (C/o — U\)n is the chemical potential. The presence of TE-TM 
splitting leads to the anisotropy of sound velocities (see [11] for details). 

(ii) U\ < 0. In this case one of the circular components is zero and the other is 
y / ne l£ '. Since there is only one broken continuous symmetry, the excitation 
spectrum consists of only one Bogoliubov branch, and the other branch is 
gaped parabolic with the gap 2|t/i|n. The chemical potential is fi = Uon 
in this domain, so that the sound velocity for the Bogoliubov excitations 
is 1 [ijrn* . 

The mean- field theory predicts an arbitrary polarization for U\ =0 since in 
this case the energy of the condensate is polarization independent. In reality, 
fluctuations destroy the order in this case at any finite temperature T . It 
can be already understood from the excitation spectrum, because, apart from 
the Bogoliubov branch, there is the gapless parabolic branch with dispersion 
h 2 k 2 /2m* and the condensate would evaporate completely due to excitation 
of these quasiparticles0 One can also map this case to the 0(4) nonlinear 
sigma model, where the order is proven to be absent for T > [T4"] . 

Note the similarity between the two-component condensates of exciton- 
polaritons and three-component condensates of spin- 1 cold atoms [121115] . Due 
to the 3D rotational symmetry, there are also only two interaction constants in 
the latter case. These constants are defined by the cross-sections of scattering 
of two atoms with the total spin and 20 Two different atomic condensates 
can also be found depending on the sign of the scattering length with the total 
spin 2: ferromagnetic and anti- ferromagnetic (or polar), which are analogs of 
circularly and linearly polarized exciton-polariton condensates, respectively. 

It is the first case, Ui > 0, that is realized in the exciton-polariton con- 
densates observed so far. The linearly polarized condensate can be seen as 



The concentration of quasiparticles with the energy e(fc) is given by 
J (2?r) -2 [exp{e(fc) /T} — l]^ 1 d 2 k and the integral diverges logarithmically for small 
k when e(k) oc k 2 . 

2 The case of the total spin 1 is irrelevant since the orbital wave function of colliding 
bosons is antisymmetric and it cannot be realized within the condensate. 



1 Mean-field description of multicomponent exciton-polariton superfluids 



5 



composition of equal numbers mutually coherent spin-up and spin-down po- 
laritons. Therefore, it is interesting to study the effect of applied magnetic field 
to this state [TS]. Considering only weak fields, when the magnetic length is 
much greater than the exciton Bohr radius, one can study only the effects of 
Zeeman field, that is described by adding 

^' = /2(|V-i| 2 -hM 2 ) (1.7) 

into Hamiltonian (11.31) . Here the Zeeman field Q is given by the half of the 
Zeeman splitting energy for a single polariton. 

To find the order parameter for the uniform condensate in the presence of 
Zeeman field it is convenient to introduce the concentrations of the compo- 
nents n±i = IV'iil 2 satisfying n+i + n_i = n. Assuming both n+i and n_i 
to be nonzero, one can take variations of the Hamiltonian 

"Hint +%' - f-m 

1 2 

= -Uo (n+i + n_i) - 2Uxn + xn-i — (fj, + f2)n+i — (fi — f2)ri-i (1.8) 
over n±i to obtain 

- 2Z7m±i = (p- U riT&)- (1.9) 
The sum and the difference of Eqs. (|1.9j) results in 

IL={U -U 1 )n, n±1 = ^(n±^j, for \fl\ < O c = n\J x . (1.10) 

For higher Zeeman fields, \Q\ > f2 c , one of the components becomes empty, 
n_i for Q > fl c and n+i for £1 < —Q c , and in this case /i = Uqu — |17|. 

Remarkably, for subcritical fields the chemical potential does not change 
at all, so that there is no change in the position of the emission line. The only 
effect of applied Zeeman field is the change of circular polarization degree g c — 
(n+i — n-i)/ (n + i +ri_i) = f2/nU\, that increases linearly with the field. The 
elliptical polarization of the condensate for subcritical fields is characterized by 
two angles, and in the the same way as for the linearly polarized condensate, 
there are two Goldstone modes; only the sound velocities change with the 
Zeeman field. This implies the full suppression of the Zeeman splitting by 
polariton-polariton interactions within the condensate [15]. Note also that 
for subcritical fields there are two phase transitions in the left and in the 
right circular component of the condensate, respectively [16] . The Zeeman 
splitting (the gap in the exciton spectrum) appears only for supercritical fields 
> fi c where the condensate becomes circularly polarized. This effect, 
observed experimentally by Larionov et al. [17], allows to measure the spin- 
dependent interaction constant U\. 
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1.3 Vortices in exciton-polariton condensates 

Vortices play a key role in various physical phenomena both on macroscopic 
and microscopic level. While the vortex formation is very important for de- 
scription of different effects in fluid mechanics, in particular, in aerodynamics 
and turbulent flow motion, the understanding of properties of quantized vor- 
tices is crucial for description of phase transitions in condensed matter. The 
well known examples are the phase transitions in type II superconductors 
in applied magnetic field, which are related to the formation and melting 
for vortex lattices [IB], and the Berezinskii-Kosterlitz-Thouless (BKT) phase 
transition P P ? l l2"T ] . 

As it was discussed above, the exciton-polariton condensates possess two- 
component order parameter (|1.2p and these condensates allow half-quantum 
vortices (half- vortices) [22]. Moreover, the half- vortices are basic topologi- 
cal excitations in this case (see [23] for a review on basic properties of half- 
quantum vortices). In spite of recent observation of both integer [51] and 
half- integer [25] vortices in exciton-polariton condensates, the presence of half- 
vortices was recently questioned [26 for the case of two-band dispersion with 
TE-TM splitting of polariton band given by Eq. (|1.4j) . In this section I present 
the details on how the vortex solutions should be found in this case (a short 
summary of this theory has been given in |27j ) . In what follows only the case 
of zero Zeeman field will be considered. 

For a 2D system of radius R the energy of a vortex is finite but logarith- 
mically large, 

E VOI =E c + E a ln(R/a), (1.11) 

where a = h/y / 2m* /j, is the characteristic radius of vortex core (the effective 
mass m* is defined below in Eq. (|1.17I) V The fact that E VOT diverges logarith- 
mically at R — > oo is good: it prevents the single vortices to be excited at 
low temperatures and thus protects the long-range order of the condensate. 
Knowledge of prefactor E s allows to estimate [21] the BKT transition tem- 
perature T c . The proliferation of single vortices appears when the free energy 
E vol - — TS crosses zero. The vortex core area is a 2 and it can be appear in 
R 2 /a 2 places, so that the entropy S = \n(R/a) 2 and this gives T c — E s /2 
if one neglects the energy of the core E c in (jl.ll|) . The energy E s \a(R/a) is 
elaborated on large distances from the vortex core r»a, which we will refer 
to as the elastic region, and the study of vortices should begin with establish- 
ing the behavior of the order parameter in this region. When TE-TM splitting 
is present this behavior is, in general, nontrivial. 

1.3.1 The order parameter on large distances 

In the elastic region the order parameter changes within the order parameter 
manyfold, i.e., the polarization of the condensate is linear everywhere in this 
domain. The circular-polarization components ij)±i defined in (|1.2p can be 
written in cylindrical coordinates (r, <fr) as 
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^ ±1 (r »a,0) = J^eW^OW (L12) 

where n is the constant concentration of the condensate at large distances, 
and the phases are written in terms of total phase angle 6 and polarization 
angle r\. These angles do not depend on the radius r (such dependence would 
only increase the vortex energy) , but they are functions of the azimuthal angle 
4>. Since the order parameters should be uniquely defined in the whole space, 
one has 

r)((t> + 2tt) -?/(» = 27rfc, 6{(f) + 2tt) - 8((f>) = 2vrm. (1.13) 

These conditions divide all possible solutions of GPE into topological sec- 
tors. Each sector is defined by two topological charges (or winding numbers), 
k and m. The state from one sector cannot be continuously transformed into 
another sector, or, in other words, any state of the condensate evolves within 
its own topological sector. The sector k — m — is the ground state sector; 
the minimum energy here is reached for position-independent order param- 
eter. By definition, the vortex is the state that minimizes the energy in a 
topological sector with at least one non-zero winding number. The energy of 
the (k, m)-vortex (|1.11[) is counted from the ground state energy, i.e., it is the 
difference between the minimal energy in the (k, m) sector and the minimal 
energy of (0,0) sector (the ground state energy). Since only the sum and the 
difference, 8 ± ij, enter Eq. (|1 . 12|) . the winding numbers can be either both 
integer or both half-integer, and the corresponding vortices are referred ac- 
cordingly. Note also that the vortex corresponds to a minimum of Hamiltonian 
H for specific boundary conditions: 5H/5xjj* = for the vortex solution and, 
therefore, it is a static solution of GPE (|1.1[) F1 

According to (|1.13|) one can add any periodic functions of <p to /?(</>) and 
#(</>) without changing the topological sector. The proper functions 7](<p) and 
0(4>) for the (fc, m)-vortex should be found from minimization of Hamiltonian 
in elastic region. The corresponding part of Hamiltonian is related solely to 
the kinetic energy term J TcPr. After substitution of (|1.12l) into (|1.4b ) and 
use of the asymptotic behavior of the complex derivative for r — > oo, 

oz 2r 0(f> 

one obtains the product of integrals over r and </> that results in the second 
term of Eq. (|1.11[) . The integral over r diverges logarithmically and should be 
cut by the core size a at small r, and by system radius R at large r. This gives 
the factor ln(i?/a). The prefactor is then given by 



3 Note, however, that this does not imply that a single vortex gives an absolute 
minimum of the H in the corresponding topological sector. For example, the 
integer vortex (1,0) can be unstable with respect to decay into the pair of (|, |) 
and (i ,— i) half-vortices for mi < m t (see subsection ll.31 2 below). 
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H 2 n f 27r 

E s = - — - / {[1 + 7 cos(2m)](1 + u') 2 + [1 - 7 cos(2u)]0' 2 }d0, (1.15) 
2m* Jo 

where the prime denotes the derivative over and 

u{4>) = r](4>) - 4>. (1.16) 

The effective mass m* and the TE-TM splitting parameter 7 are defined in 
(CI3 by 

11/1 1 \ m t — mi . , 
= 77 — + — , 1=— (1.17) 

m* 2 \mi m t J m t + mi 

Variations of the functional (I1.15[) over and u lead to the equations 

[1 - 7Cos(2w)] 0" + 2-/sm(2u)u'9' = 0, (1.18a) 

[l + 7COs(2M)]u" + 7 sin(2u)(l-M /2 -6l' 2 ) =0. (1.18b) 

In general, the polarization will be radial at least at one specific direction 
and it is convenient to count the azimuthal angle from this direction and set 
the total phase to be zero at this direction as well. Then, the solutions of 
Eqs. (jl.lSa .b 1 ) for (k, m)-vortex should satisfy the boundary conditions 

u(Q)=Q, 0(0) = 0, (1.19a) 

u(2tt) = 2(k - 1)tt, 0(27r) = 2m7r. (1.19b) 
The solutions in question are trivial for some particular vortices. 

(i) Hedgehog vortices. These are (1, m)-vortices having = m<p and u = 0, so 
that the polarization angle 77 = cf>. Polarization points into the radial di- 
rection everywhere and these vortices look like hedgehogs. These solutions 
are similar to magnetic monopoles |28j . 

(ii) Double- quantized polarization vortex (2,0). In this special case = 0, 
but u = (f>, resulting in r\ = 2cj>. Polarization rotates twice when one 
encircles the vortex core. These vortices and experimental possibilities of 
their excitation in exciton-polaritons fields have been studied by Liew at 
el. Ell- 



in other cases the solutions should be found numerically. Both Eqs. (|1.18k ) 
and (|1.18b ) can be integrated once to give 

[l-7cos(2tO]0' = Ci, (1.20a) 

[1+7 cos(2u)] u' 2 + [1 - 7 cos(2u)] 0' 2 - 7 cos(2u) = C* 2 , (1.20b) 

and the solutions can be written as integrals of elementary functions. The 
constants should then be found, e.g., by shooting, to satisfy the bound- 
ary conditions (|1.19l) . The functions 0(</>) and rj((f)) are shown in Fig. 11.11 for 
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Fig. 1.1. The dependence of polarization angle r\ (solid lines) and phase angle 6 
(dashed lines) on the azimuthal angle <f) for two values of TE-TM splitting parameter: 
7 = —0.4 (thin lines) and 7 = —0.9 (thick lines). The panels show the behavior of 
the angles for the (3,3) half-vortex (a), the ( — 313) half-vortex (b), the ( — 1,0) 
polarization vortex (c), and the (0, 1) phase vortex (d). In the last case the periodic 
function rj{4>) has been upscaled for clarity 



elementary half-vortices and for two integer vortices (—1,0) and (0,1), that 
also exhibit nonlinear dependencies of polarization and phase angles. 

Fig. 11.11 demonstrates the behavior of angles for negative values of the 
TE-TM splitting parameter 7. The functions 9(<f>) and ?y(</>) for positive 7 can 
be found by the shift. Indeed, the change u — » u + (tt/2) in Eqs. (|1.18b .b) 
results in the change of the sign of 7. More precisely, to satisfy the boundary 
conditions (JTTT9J) the transformations can be written as 



7 -7. 

u{4>) ->• u 1 4> + 



6(4,) 



2{k 

IT 



TV 

2' 



2(k-l) 



2(k- 1) 



(1.21a) 
(1.21b) 

(1.21c) 



and they can be applied to all vortices except the hedgehogs with k — 1 (and 
where they are not necessary, of course, since u(4>) = 0). 
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The nonlinear change of angles seen in Fig. II .11 becomes especially evident 
when 7 approaches ±1. This limit correspond to a strong inequality between 
effective masses, e.g., m; 3> m t for 7 — > — 1. Qualitatively the strong nonlin- 
earities can be understood if one introduces the effective masses for the phase 
me and for the polarization m^, 



1 cos 2 m sin 2 u 1 sin 2 u cos 2 u 



mg m t mi m v m t mi 

and writes the energy (|1.15p as 



ti 2 



n 

T 













m e J 



(1.22) 



(1.23) 



The effective masses (|1.22l) depend on the orientation of polarization. Since 
u(</>) changes between the values specified by (|1.19p . there are sectors where 
me ~ mi and m v sa m t , and there are sectors where me « m t and m v « mi . 
To minimize the energy (| 1 . 23[) in the case mi ^ m t , the phase angle changes 
rapidly and the polarization angle stays approximately constant in the former, 
while there is the opposite behavior in the latter. 

1.3.2 The energies and interactions of vortices 

In the absence of TE-TM splitting the energies E s of vortices are 

E§ m) =E (k 2 +m 2 ), E = ^, (for 7-O). (1.24) 

It is seen that in this case the energy of an elementary half- vortex is exactly 
half of the energy of an elementary integer vortex. Important consequences 
can be drawn from this relation concerning the interactions between half- 
vortices. The four elementary half-vortices can be divided in two kinds, right 
half- vortices with k + m = ±1, and left ones with k — m = ±1. One can see 
from (|1.12[) the the right half- vortices possess the vorticity of the left-circular 
component of the order parameter, the amplitude of this component goes to 
zero and the phase of this component becomes singular in the vortex core 
center, and, as a result, the polarization becomes right-circular in the core 
center. For left half-vortices the picture is opposite. It follows from (|1.24D 
that the left and the right half-vortices do not interact with each other. Con- 
sider, for example, the (5,5) and (—5,5) half- vortices. The elastic energy 
of this pair is Eohi(R/a) both when they are far away from each other and 
when they are in the same place forming the phase vortex (0, 1). So, the only 
possible coupling between the left and the right half- vortices is of short range, 
related to the overlap of their cores and resulting change of the core energy 
term E c in (jl.lip . So, the only long-range coupling is present between the 
half-vortices of the same kind. It can be shown that identical half-vortices 
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repel each other logarithmically, while the half-vortices and anti-half- vortices, 
(k,m) and (—k,—m), attract each other logarithmically, as it is in the case 
of vortices and antivortices in one-component condensates |30j . This simple 
picture is changed in the presence of TE-TM splitting that leads to the long- 
range interaction between half-vortices of different kind. 

The logarithmic prefactors Es k for elementary half-vortices and ele- 
mentary integer vortices are shown in Fig. 11.21 For all of them, expect the 
hedgehog (l,0)-vortex, these energies are even functions of 7, which can be 
proven using the transformations (| 1 . 2 1 a|) . These energies decrease with in- 
creasing 7 2 . The case of hedgehog is special, as it has been discussed above. 
The hedgehog polarization is radial everywhere and Eg 1 is defined purely 
by the longitudinal effective mass, 

E^ 1 ' ^ = Trh 2 n/mi = E (l +7). (1.25) 

It can be shown that when two vortices with winding numbers (fci,mi) and 
(A;2, m-i) are injected in the condensate and they are separated by distance r, 
such that a <C r <C R, the energy of the condensate is increased in logarithmic 
approximation (i.e., omitting the core energies) by 

^(fci+fca.mi+ma) ln (fl/ a ) 

+ \ E (k umi ) + E (k 2 ,m 2 ) _ E (k 1 +k 2 ,m 1 +m 2 )l m ( r / fl )_ ( L26 ) 

The second term in (jl.26[) gives the interaction energy of two vortices. The 
coupling between vortices arising due to TE-TM splitting can be analyzed 
analytically in the limit of small 7. 

The solutions of Eqs. (|1.18a .b) for k ^ 1 are written as series in 7, 

7sin[2(fc - 1)4>] 



2(jfe-l 



m\{k — l) 2 — (1 — m 
16(fc- l) 3 



•2\ 



7 2 sin[4(fc- 1)0] + (1.27a) 



»(.0 = ( /,-l) f) - t( fc V* k + J-™1 1)] 7si„[2(/.- Ho] 



+ [5(fc - 1)4 + 2( ^f (TO ;- 3) + (m2 ' 1)2] 7 2 sinWfc - m 
64(k — l) 4 



Substitution of these expression into (j!.15[) gives 



(1.27b) 



E^ m) _ 2 2 [fc 2 (fc-2) 2 + 2[2 + 3fc(fc-2)]m 2 +m 4 ] 2 _ 
Eo ~ 1 ' 8(fc-l) 2 7 1 } 
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TE-TM splitting parameter 7 

Fig. 1.2. The logarithmic prefactor of vortex energies E s (see Eqs. 1)1.110 and (|1 . 15|) ) 
for half-vortices and integer vortices as functions of TE-TM splitting parameter 7 
Q1.170 . The curves are labeled by the winding numbers (k,m) of the vortices, and 
the energies are given in the units of Eg = irh 2 n/rn* 



and, in particular, 



Ei- 1 ® = En 



1 7 

32 ' 



El * 1 ) = E 



(1.29) 



There is no difference between the energies of half- vortices at this order of 
7-serics. The difference, however, appears in the next order. The series for the 
angles up to 7 4 are rather cumbersome to be presented, but they result in 



E, 



E' 



2' ± 2, 



9 1 9 . 



Eo 
2 

Eo 
2 



2 16 



r_ 
2 



117 4 

144 



(1.30a) 



(1.30b) 



Eqs. (|1.25|) . (|1.29l) . and (jl.30p can be used to find the interactions between 
half- vortices according to Eq. (|1.26|) . Most important interaction that appears 
due to TE-TM splitting is between the (|, 5) and (i, —5) half- vortices. For 
small 7 their coupling constant is linear in 7 and the interaction energy is |31) 



-) ('- 



-jE \n(r/a). 



(1.31) 



It should be noted that the interrelation between m t and to/, i.e., the sign of 
7, depends on the detuning of the frequency of the cavity photon mode from 
the center of the stop-band of the distributed Bragg mirror [32]. So, one can 
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have both attraction and repulsion of the (5,5) and (I, — i) half- vortices. The 
coupling of the other left and right half-vortices is quadratic in 7. The (— |, i) 
and (— i) — |) half- vortices repel each other with the interaction energy being 
-(7/32) 7 2 £ , ln(r/a). The (-|,±±) and (|,±|) half-vortices attract each 
other with the interaction energy being {1/%)^ 2 E ln(r/a). 

In the absence of TE-TM splitting there is no coupling between the right 
half-vortices (with km > 0) and the left ones (with km < 0) and there are 
two decoupled BKT transitions, corresponding to the dissociation of pairs of 
left and right half- vortices [2HIS]- The transition temperature is then esti- 
mated from the energy of single half- vortex as Eq/A. The TE-TM splitting 
of polariton bands changes this picture substantially. First, because all four 
half-vortices become coupled and, secondly, because the energies of a vortex 
and its antivortex become different, so it is not clear with one should be used 
in the estimation of critical temperature. 

One expects qualitative modifications of the BKT transition in the region 
of 7 close to —1. In this region the attraction of the (5,5) and (5,-5) half- 
vortices becomes very strong and, as a result, the hedgehog is the vortex with 
the smallest energy in the system for 7 < j c ~ —0.6 (see Fig. II. 2p . It does 
not mean, however, that the transition temperature can be estimated from 
the energy of the hedgehog in this region. In fact, the phase transition occurs 
due to dissociation of vortex-antivortex pairs, and the energy of the (1, 0) and 
(—1,0) pair is still bigger than the energy of the pair of two half- vortices. 
One expects that when pairs of half-vortices are thermally excited in the 
system they will tend to form molecules consisting of the hedgehog (formed by 
merging of the (i, \) and (|, — |) half-vortices) with the (— |, — |) and (— |, i) 
half- vortices being attached to it. The proliferation of these (— i, — |) — (1, 0) — 
(—5, 5) molecules defines the phase transition for mi 3> m t . 



1.4 Geometry of the half- vortex fields 

In general, two coupled Gross-Pitaevskii equations for the circular components 
of the order parameter 

+ [U (\iJj +1 \ 2 + |V;-i| 2 ) - 2(7 1 |^_ 1 | 2 ] (1.32a) 



+ [£/o(|^-i| 2 + |^ + i| 2 ) - 2[/ 1 |^ +1 | 2 ] ^_ X) (1.32b) 

are not separated in cylindrical coordinates (r, </>) . The variables are sepa- 
rated only in special cases of hedgehog vortices and the double-quantized 
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Fig. 1.3. Showing the geometry of half-vortices for different values of TE-TM split- 
ting parameter 7 (JTTT7J). The plots are obtained from numerical solutions of GPEs 
(|1.32J1 [33] . The interaction constants are related as Ui = 0.55C/o- The local polar- 
ization ellipses are drawn with the thick (red) lines. The streamlines of the current 
are shown by thin (green) lines. The panels demonstrate the following cases: (a) the 
half- vortices (~,±i) for 7 = —0.5 (the lemon morphology); (b) the half- vortices 
(|,±|) for 7 = 0.5 (the lemon morphology); (c) the half-vortices (— i,±i) for 



7 = —0.5 (the star morphology); (d) the half-vortices (- 
morphology) . 



, ±|) for 7 = 0.5 (the star 



polarization vortex discussed in previous section after Eqs. (jl,19[) . For other 
vortices one needs to solve GPEs numerically in two spacial dimensions with 
the boundary conditions at large distances defined by Eqs. (|1 . 12|) and (jl.181) . 

Each vortex with nonzero phase winding number m is characterized by a 
finite superfluid current circulating around its core center. Performing numer- 
ical solutions it is important to take into account the fact that streamlines 
of the current are deformed with respect to perfect circles in the presence of 
TE-TM splitting. The circular components of the current J are given by 
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^(^-*--.^)}. <-> 

and J_i = Jii- They are related to the radial J r and the azimuthal 
components of the current by 

J r = (e i0 J+i + e-**J_i) , J = -±= (e i0 J + i - e^J-i) . (1.34) 

At large distances r a — h/ ^/2m*/x one can use Eq. (|1 . 12|) to find 

hn dO 
J r = — — 7sm(2u) — , 1.35a 

J* = — [1-7«m(2«)]^. (1.35b) 

Note that the condition of conservation of the total number of polaritons for 
the static vortex solution of Eqs. (|1.32[) . 



div J = - 

r 



0, (1.36) 



implies dJ^/dcj) = 0. So, Eq. (|1.18b ) obtained in the previous section is in fact 
the condition of conservation of the azimuthal current. 

The warping of streamlines of current is shown in Fig. 11.31 The order 
parameter has been found numerically 33 for different values of TE-TM 
splitting parameter. To find the static solutions of GPE (|1.32j) . we have been 
choosing an initial order parameter ip(r,t = 0) satisfying the boundary con- 
ditions that follow from Eqs. (|1.18[) and that are shown in Fig. 1 1.1 1 for a given 
topological sector (k,m). Apart from this the initial functions were contin- 
uous but arbitrary. Then the functions were evolved according to (|1 .32|) in 
imaginary time. As a result, the order parameter relaxed to corresponding 
static half- vortex solution. The resulting half-vortex solutions are found to be 
independent of the initial shape of ip(r, t = 0). 

In Fig. 11.31 one can see two distinct morphologies of basic half- vortices. 
The geometry of half-vortex solutions can be discussed in terms of singular 
optics 34 , 35j , where the polarization singularity related to a half- vortex is 
referred as C-point, to indicate that the polarization is circular at the vortex 
center and, therefore, the direction of the main axis of polarization ellipse 
is not defined. The morphologies of the field around C-points are classified 
by the index of associated real tensor field, and, additionally, by the number 
of strait polarization lineal that terminate at C-point |36j . The tensor index 



4 The tangents of polarization lines define by the direction of the main axis of 
polarization ellipse in each point. 
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coincides with the polarization winding number k, and the number of lines 
could be either one or three. As a result, three different morphologies can be 
found. Following Berry and Hannay [37J[35], these morphologies are referred 
as lemon, star, and monstar. 

The lemon configuration is characterized by k = | and by only one straight 
polarization line terminating in the vortex center. This is the morphology of 
vortices in Fig. I1.3h ,b with the straight polarization line being defined by 
= 0. The star configuration is characterized by k = — i. In this case there 
are always three straight lines terminating in the vortex center. The stars 
are realized in Fig. I1.3b .d and three straight polarization lines are defined by 
(/) = 0,±2tt/3. 

The change of parameter 7 leads to deformation of polarization texture 
and to deformation of streamlines of the current, but it does not result in 
the change of morphologies of half- vortices. In principle, one could expect 
the transformation of lemon into monstar, since these morphologies possess 
the same topological index k — i. Contrary to the lemon case, however, the 
monstar is characterized by three straight polarization lines terminating in the 
vortex center, similar to the star configuration. So, from geometrical point of 
view the monstar has got intermediate structure between the lemon and the 
star, and this is why its name is constructed from "(le)mon-star" . 

To have the monstar configuration one needs a special behavior of polar- 
ization angle ??(</>). Namely, it is necessary to have 



drj 



> 1 for k = -. (1.37) 

4>=o ^ 



In this case the polarization angle initially rotates faster than the azimuthal 
angle <p, but since the total rotation of 77 should be still 7r when <j) is changing 
up to 2tt, as it is dictated by the winding number k = i, there will be three 
roots of the equation r]((f>) — <f>. These roots, and ±^ m , define three straight 
polarization lines terminating in the vortex center for the monstar geometryQ 
One can see from Fig. ll.la and ll.3a that when 7 approaches —1 the derivative 
becomes very close to unity, but it never becomes bigger than 1, so that the 
monstar is not formed. The reason preventing the appearance of the monstar 
is that it is not energetically favorable to satisfy the condition (jl.37[) . In fact, 
it is the most energetically favorable to have r\ = <fi, as for the hedgehog — the 
vortex having the smallest energy when 7 — > — 1 (see Fig. ll.2|l . The rotation of 
polarization of the half- vortex is also synchronous with the azimuth in rather 
wide sector, but the polarization never overruns the azimuth. The monstar 
half- vortices, however, are expected to be found in the exciton-polariton con- 
densates out of equilibrium [39,9,40 , where their appearance is not restricted 
by energetics. 



5 Note that for the monstar all polarization lines residing within the sector — <f> m < 
<f) < 4>m terminate in the vortex center, but only three of them are straight, i.e., 
are having nonzero inclination at r — > (see [36] for the details). 
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1.5 Four-component exciton condensates 

Excitons formed by an electron and a heavy hole in the semiconductor quan- 
tum wells can be in four spin states [41]. The states with the total spin pro- 
jection ±1 are optically active. These bright excitons are formed by the heavy 
hole with the spin +| and the the electron with the spin — i, or by the — f 
hole and the +^ electron. The other two states are hidden from the observer 
and are usually referred to as the dark excitons. The total momentum of these 
states is ±2 and they are formed either by the +§ hole and the +\ electron, 

Q 1 Z Z 

or by the — § hole and the — ^ electron. 

The exciton-polaritons discussed in the previous sections are coupled states 
of quantum-well excitons and microcavity photons. Only the bright exci- 
tons are involved in this coupling, and the resulting condensates are two- 
component. Since the frequency of a single exciton-polariton is shifted down 
with respect to the single exciton frequency by a half of the Rabi frequency, the 
presence of dark excitons is irrelevant in this case provided the exciton-photon 
coupling is strong enough. Contrary, when pure exciton condensates are of in- 
terest, all four exciton spin states should be, in general, taken into account. 
The formation of exciton condensates is possible for cold indirect excitons 
in coupled quantum wells [321 [3211111113 US]. The life-time of these excitons 
is long enough, the excitons can travel coherently over long distances, and 
the condensates can be formed in quasi-equilibrium conditions. The presence 
of four-component exciton condensates has also been experimentally demon- 
strated recently [47] . 

The indirect excitons are dipoles oriented along the growth axis of the semi- 
conductor structure, and their main interaction is spin-independent dipole- 
dipole repulsion. The condensate state, however, is defined by weak spin- 
dependent interactions arising from electron-electron, hole-hole, and exciton- 
exciton exchanges [48]. In what follows, I will assume the signs of these inter- 
actions to be such that they favor the distribution of excitons over all four spin 
states, populating both bright and dark components. This state is similar to 
the linearly polarized two-component condensates described above, but there 
is one important qualitative difference between them. The exchange scattering 
of two excitons can result in transformation of their spin states [49] . Namely, 
two bright excitons can turn into two dark ones after collision and vice versa. 
These processes are described microscopically by the Hamiltonian 



Umix — W 



V'+2^L 2 '0+iV'-i + ^+1^-1^+2^-2 ■ (1-38) 



In mean-field approximation, the creation ip^ and annihilation ij) a exciton 
operators (a = ±1, ±2) are replaced by the order parameter components, ip* 
and ?/v, respectively. The contribution of the resulting exciton- mixing term 
Hmix into the total energy of the exciton condensate depends on the rela- 
tive phases of the components. The term of this type is absent in the two- 
component exciton-polariton case. Remarkably, the mixing of excitons always 
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leads to the decrease of the condensate energy, which is achieved by hxing the 
proper interrelation between the phases. Denoting by 9 a the phase of ifj a , one 
can see that the following relation holds within the order parameter manifold 



The mixing term (|1.38[) additionally favors the formation of the four- 
component exciton condensate with equal occupations of the components. 
This fact can be seen from different perspective. H m i x describes the transfor- 
mation of pairs of excitons, and, in the same way as in the BCS theory of 
the superconductivity, this term leads to the pairing of particles. This pairing 
leads to a decrease in the energy of the system and results in appearance of 
the gap in the excitation spectrum for one excitation branch. The other three 
excitation branches are Bogoliubov-like. This follows from the fact that the 
phase locking condition (|1.39|) leaves there angles to be undefined, so that 
there are three Goldstone modes apart from the gaped mode induced by the 
mixing. 

The effect of applied Zeeman field on four-component exciton condensate 
is expected to be very spectacular [35]. The Zeeman splitting is different for 
dark and bright excitons: the g-factor is given by the sum of the electron and 
hole (/-factors for the former, and by their difference for the latter. When such 
a field is applied to the exciton condensate its action is two-fold. On the one 
hand, it polarizes the bright and dark components with different degrees of 
circular polarization, and thus reduces the Zeeman energy of the condensate. 
On the other hand, the induced imbalance in the occupation of the components 
increases the energy of the mixing term (|1 .38[) and suppresses the gap in the 
spectrum discussed above. The interplay between these two effects can lead to 
the first-order transitions from the four-component exciton condensate to the 
two- or the one-component condensates. Note also that due to the presence of 
'Hmix the system of equations defining the concentrations of the components 
of the exciton condensate in the Zeeman field is nonlinear, contrary to the 
case of two component exciton-polariton condensate (see Eq. (|1.9JI ). 

Finally, it is important to note that vortices in the four-component exciton 
condensate in the presence of the mixing of the component are composite: the 
vorticity of one component should be accompanied by the vorticity of another 
component to satisfy the the phase-locking condition (jl.39l) . As a result one 
expect twelve elementary vortices. These are four polarization vortices (two 
in each components), and the eight paired half- vortices in the bright and dark 
components. 

1.6 Conclusions and perspectives 

The mean-field approximation provides simple and reliable method to study 
the polarization properties and topological excitations of exciton-polariton 




(mod2vr), if W < 0, 
(mod 2tt), if W > 0. 



(1.39) 



1 Mean-field description of multicomponent exciton-polariton superfluids 



19 



and exciton condensates that possess two and four components of the order 
parameter, respectively. This includes, in particular, the description of the 
polarization of the ground state and elementary excitations of the condensates 
and their change in applied Zeeman field, as well as the description of the 
texture of vortices and vortex interactions. 

The elementary topological excitations in two-component exciton-polariton 
condensates arc four half- vortices (k,m) with k,m — ±|, characterized by 
half-quantum changes of polarization and phase angles. In the absence of 
transverse-electric-transverse-magnetic (TE-TM) splitting of the lower polari- 
ton band there is no coupling between the left half-vortices (with km < 0) 
and the right ones (with km > 0), and one expects two decoupled Bcrezinskii- 
Kosterlitz-Thouless (BKT) supcrfluid transitions happening at the same tem- 
perature in the system. The TE-TM splitting results in two qualitative effects. 
First, the cylindrical symmetry of the half- vortex field is spontaneously broken 
that leads to warping of the polarization field around a half- vortex and to devi- 
ation of the streamlines of the supercurrent from the perfect circles. Secondly, 
there appears long-range interactions between left and right half-vortices. 
These interactions are particularly important in the case of large longitu- 
dinal polariton mass mi when it favors the formation of hedgehog (monopole) 
vortices (1,0) from the (|, |) and (|,— \) half-vortices. The peculiarities of 
the supcrfluid transition in this case and related features of polarization tex- 
tures of the exciton-polariton condensates are subjects of further studies. In 
what concerns the geometry of the half- vortex field it is shown that only two 
configurations, lemon and star, are realized. The monstar configuration is not 
energetically favorable for any value and sign of TE-TM splitting. 

The essential feature of four-component exciton condensates is the pres- 
ence of mixing and related phase locking between dark and bright excitons. 
One expects a nontrivial Zeeman-field effect resulting in a discontinuous 
change of the polarization state of the condensate in course of the first-order 
transition. The presence of composite vortices in different components should 
lead to the formation of interesting polarization patterns in driven exciton 
condensates that provide an important topic for investigation, both experi- 
mental and theoretical. 
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